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ABSTRACT 

We investigate the X-ray properties of three interacting luminous infrared galaxy systems. 
In one of these systems, IRAS 18329+5950, we resolve two separate sources. A second, IRAS 
20550+1656, and third, IRAS 19354+4559, have only a single X-ray source detected. We 
compare the observed emission to PSF profiles and determine that three are extended in emis- 
sion. One is compact, which is suggestive of an AGN, although all of our profiles have large 
uncertainties. We then model the spectra to determine soft (0.5-2 keV) and hard (2-10 keV) 
luminosities for the resolved sources and then compare these to relationships found in the liter- 
ature between infrared and X-ray luminosities for starburst galaxies. We obtain luminosities of 
log(L so f t /L ) = 7.32, 7.06, 7.68 and log(L ha rd/L ) = 7.33, 7.07, 7.88 for IRAS 18329+5950, 
IRAS 19354+4559, and IRAS 20550+1656, respectively. These are intermediate to two sep- 
arate predictions in the literature for star-formation-dominated sources. Our highest quality 
spectrum of IRAS 20550+1656 suggests super-solar abundance of alpha elements at 2a signif- 
icance, with log = M - 0-4 ± 0.2. This is suggestive of recent enrichment with Type II 
supernovae, consistent with a starburst environment. The X-ray properties of the target galaxies 
are most likely due to starbursts, but we cannot conclusively rule out AGN. 
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1. Introduction 



Over the past few decades, much work has emphasized the role and importance of mergers in shaping 
galactic evolution ( Toomre & Toomre 1972 : Ellison et al. 2011 : Hopkins et all2006 ). In particular, mergers 



have been suggested as potential drivers of starbursts and AGN activity (Bauer et al.ll2004l : ISummers et al. 



20041 : iBell et al ll201ll : lKomossa et alJl2003l:lTeng et al.l2009l : lpiconcelli et alhoiol) . Broadly, the hierarchi- 
cal growth simulations in lHopkins et all (12006 ) paint a picture in which galactic mergers provide gas inflows 
to ignite starbursts, which are then followed by growth of the central supermassive black holes (SMBHs) of 
galaxies, sometimes leading to the birth of an AGN until feedback possibly stops the process. The interac- 
tions of two gal axies can also c ause both of their central SMBHs to become active simultaneously, resulting 
in a dual AGN (|Komos sa 2006). A complication in this scenario is that some simulations suggest that most 
of the dual AGN activity is non- simultaneous, particularly dur ing their highest luminosity phases when they 
would also be easiest to detect (|Van Wassenhove et al.ll2012r) . Other simulations have suggested that this 
merger scenario was more important in the past, when interacting sys tems were much more gas -rich, but 
that it still plays a role in the local universe for lower-luminosity AGN ([Draper & Ballantyndl2012l) . 



Lumin ous infrared galaxies (LI RGs), galaxies that have infrared luminosities Ljr > 10 L (for a 



Sanders & Mirabel 1 19961 and references therein), where IR ranges here span 8-1000 /im, are 



review, see 

often found to be mergers and have distinctive morphologies (|Yuan et al.ll2010l) . The emission from such 
systems is prominent in the infrared (IR) portion of the spectrum due to heavy dust reprocessing. The Great 
Observato ry All-sky LIRG S urv ey (GOALS), which is observing LIRGs acro ss the electromagnetic spec 
trum (e.g. 



Armus et al 



20091 and 



2009; lAlonso-Herrero et al 



Iwasawa et al 



2011 



201lb. and o ther studies (e.g. IPtak et al.ll2003 



Charmandaris et al. 



Braito et al. 



2010l) have attempted to determine if the drivers of 



the high IR luminosities in these systems are central starb ursts or AGN. In one s uch system, NGC 6240, two 
separate AGN from merging galaxies, were detected by iKomossa et al. ( 2003 ). This discovery has fueled 
much discussion about the nature of dual AGN, their prevalence, and lifetimes. 

Several st udies have investigat ed the connection between AGN activity and environment. When look- 
ing at quasars, ISerber et al. 1 J2006h found a local overdensity of ~3 around the brightest objects and an 
overdensity of ~ 1.4 for fainter quasars. They also found that t he overdensity aroun d quasars relative to L* 
galaxies is about unity at 1 Mpc, but increases at smaller scales. Ellison et al.l (1201 lh reported a similar result 
focusing on AGN; they found AGN are 2.5 times more likely to be in pairs than sim ilar but inactive galax- 
ies. Using the Swift's Burst Alert Telescope (BAT) AGN sample, iKoss et al.l (1201 lh also found evidence 
that AGN activity is associated with mergers and interactions. Specifically, they report that 24% of BAT 
AGN are undergoing a merger, whereas this number is closer to 1% for a sample of normal galaxies. These 
works and work done by other authors (e.g.. lAjeilo et al 1 l2012l) suggest that a galaxy's environment can play 
a non-negligible role in fueling AGN, whether through interactions or mergers with nearby galaxies. 
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Because of heavy obscuration at most wavelengths in these dust-enshrouded systems, X-rays, which 
are less attenuated, become a useful tool to probe behind optically thick screens. The hard X-ray band (2- 
10 keV) luminosity of AGN is typically described by a power law with photon index T ~ 2. A proposed 
explanation for this spectral behavior is that an accretion disk feeding the central SMBH thermally emits 
in the opt ical and UV. These photons inverse Compton upscatter to X-rays off of a h ot plasma surrounding 
the disk (|Shapiro et all 1 19761 : Isunyaev & Titarchuklll98ol : Haardt & Maraschilll993l). Much of a galaxy ' s 
SMBH growth is believed to occur in regions heavily obscured by dust (|Hopkins et al.ll2006l : lFabianl ll999). 
so searching in X-rays is one of the few viable options for spotting this phase, provided the column density 
is not overly large (> 10 24 cm -2 ). 

AGN activity is not the only interesting phenomenon strongly associated with galactic interactions. 
Another prominent event often triggered is a marked increase in the star formation rate of one or both 
compone nt galaxies, provided neither is gas poor. If this increase is large enough, the galaxy becomes a 
starburst. Muzzin et all (|2012h suggest that, while the gross properties of star-forming galaxies are tied to 
their stellar masses, their environment appears to regulate the fraction of systems that are starbursting. 



Like AGN, star-forming regions also produce specific X-ray signatures (e.g. jGeorgakakis et al.ll201 lh - 
The ionizing photons from supernovae and short-lived O and B stars form a thermal plasma typically be- 
tween 0.1-1 keV. A secondary po wer law compon ent, associated with high mass X-ray binaries (HMXBs), 
indicative of recent star formation (|Fabbianoll2006r) . is often observed in starburst systems. These are said to 
trace recent or ongoing star formation because they require an OB star around a neutron star (NS) or black 
hole (BH) and thus only have a lifetime of 10 6-7 years. HMXBs have been observed with a range of photon 
indices, typically lying within 1-1.5, while low mass X-ray binaries (LMXBs), a possible contaminant, are 
usually described by ei ther a power law with index ~1.6 or bremsstrahlung at 7.3 keV (IFabbiand 120061 : 
Persic & Rephaelill2003T) . These spectral indices can run as steep as 2-3 in some cases. LMXBs do not trace 
recent star formation, however, since they are associated with lower mass companions around a NS or BH 
and thus have lifetimes that are 100-1000 times longer. 



Kennicuttl (|1998T) associated star forming regions with IR emission due to dust reprocessing. iRanalli et al 
((2003J), hereafter R03, proposed an X-ray/IR luminosity relationship in IR-bright galaxies, arguing that star 
formation is traced in the IR from reprocessing and in X-rays from HMXBs, provided there is no AGN 
present. The use of X-rays as a tracer of star formation also requires that LMXBs are a negligible con- 
tributor, which is generally a safe assumption for the high star formation rates encountered in starbursting 
systems. One can then use these relations to predict the X-ray output of a system originating from star 
formation alone . It should be noted, however, that th e existence of either a starburst or an AGN does not dis- 
count the other (ISani et al .1120 1 (j : I S antini et al.ll2012h . Additionally, some argue for a more complete picture 
taking int o account escaping, n on-reprocessed UV radiation when attempting to synthesize star formation 
data (e.g Jvattakunnel et al.lEoilh . 



In this paper we attempt to use the association between AGN activity, starbursts, mergers, and, by 
extension, large IR luminosity, to search for dual AGN. We also study the X-ray/IR relationships to better 
understand the current star formation rate due to interaction processes in LIRGs. In ^2j we describe our 
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sample selection and discuss the systems we examined in detail, ^presents image and spectral analysis as 
well as the luminosities we obtain. $4] compares our data to existing IR/X-ray relations and discusses our 
prospectives on dual AGN detection. Finally, we conclude in §[5] 

We adopt a cosmology of Hq = 70 km s -1 Mpc -1 , $7m = 0.27, and f^A = 0.73. 



2. Sample 

We derived our sample from a subset of nearby IR bright galaxies presented in lArribas et all ((2004). 
Our systems were nearby (redshifts of 0.03 < z < 0.07) to ensure a higher number of photon counts and 
were known to be interacting with two optical nuclear regions, though there was some uncertainty about the 
location of the secondary nucleus in IRAS 20550+1656. All were LIRGs, and many had roughly meas ured 
nuclear s eparations. We selecte d based on interaction class using the scheme proposed in lsurace ( 1998 ) and 
refined in lVeilleux et al.l(|2002r ). restricting our sample to class Ilia and Illb systems with nuclear separations 



around 10 kpc. A class Ilia system is composed of two or more galaxies with evidence of tidal interactions 
(tails, bridges) and with nuclear separation of greater than 10 kpc, whereas a class Illb system has similar 
tidal structures but with nuclear separation less than 10 kpc. Additionally, we required that systems have 
angular separations that would likely be resolvable with XMM Newton, or ~ 8-10". After this, we selected 
the LIRGs with the highest IR fluxes, and, by association, likely the most X-ray counts barring the system 
being Compton thick. Three systems met our criteria: IRAS 18329+5950, IRAS 19354+4559, and IRAS 
20550+1656. We observed these with XMM Newton and restricted all analyses to 0.3-10 keV so as to stay 
in range of the detector's optimal sensitivity. 



Our observations are summarized in Table [TJ 



2.1. IRAS 18329+5950 



IRAS 18329+ 5950 is a class Ilia system composed of two merging galaxies with estimated nuclear 
separations of 28" (|Arribas et al.ll2004r ). The offset in declination (5. "70) of the two galaxies in this system 
is much smaller than the offset in right ascension (27 "5). In projection, it appears that the nucleus of the 
eastern component (18329E) is close to the plane of the disc of the western component (1 8329W). The 



system is classified as a LIRG in the IRAS Revised Bright Galaxy Sample (hereafter RBGS uSanders et al. 



2003). Using the RBGS fluxes, we calculate luminosities of log(LpiR/L ) = 11.48 (where FIR denotes 
40-400/um) a nd log(LTR/Lp) = 11.60 based on the redshift of z=0.029, which corresponds to a distance of 
~ 128.5 Mpc ( Armus et al.ll2009l) . Optical line ratios, more extensively discussed in 94.51 suggest that this 
system is a starburst or starburst/AGN composite. 
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2.2. IRAS 19354+4559 



IRAS 19354+4559 is a class Illb merger system with nuclear separations of approximately 8." 5 (lArribas et al 



2004). In this case, the offsets in right ascension and declination are comparable, being 6."7 and 5."3, respec- 
tively. Both galaxies have significant tidal features suggestive of mutual interaction and appear as edge-on 
disks, although the eastern component has several foreground stars contaminating the optical images, par- 
tially confusing the orientation. Using fluxes from the RBGS, we calculate a FIR luminosity log(LpiR/Lo) 
= 11.06 and an upper limit IR luminosity of logfLyR /L ffi ) = 1 1.8 5 based on the redshift of z=0.067, which 
corresponds to a distance of ~289 Mpc (|Lawrence et al.lll999l) . We did not find any IR or optical line 
diagnostics for this system. 



2.3. IRAS 20550+1656 



Similarly, IRAS 20550+1656 is a class Ill b merger system. The separation between the optical nu- 
clei is uncertain, but believed to be about 11: "1 (|Arribas et al.ll2004r) . From the RBGS fluxes, we calculate 
luminosities of log(LFTR / L ffi ) = 11.66 and Ios CLtr/L^) = 11.94 for a reds hift of z =0.036, correspond ing to 
a distance of ~ 161 Mpc (|Armus et al.ll2009h . Optical (|Baan et all 1 1998b and IR (|lnami et allboioh spec- 
troscopy and line ratio diagnostics indicate that this system is a starburst. The Spitzer data in particular 
(Inam i et al suggests that 80% of the IR luminosity comes from a region outside the nuclei of the two 

interacting galaxies in this merger. 



3. XMM-Newton Observations 
3.1. Image Analysis 

All data were taken with the EPIC instruments on board XMM Newton and reduced using XMM SAS 
Version 12.0.1. As is typical in X-ray image analysis, we filtered the event files to discard any time interval 



Table 1 . Observational Data 



IRAS Name 


Observation ID 


Date 


Exposure 


Distance 


Separation 


Separation 


HR 


Source Counts 








(ks) 


(Mpc) 




(kpc) 






18329+5950E 


0670140301 


2011-05-04 


9.4 


128.5 


28" 


17.4 


-0.52 ±0.20 


170 


18329+5950W 


0670140301 


2011-05-04 


9.4 


128.5 


28" 


17.4 


-0.49 ±0.1 8 


173 


19354+4559 


0670140501 


2011-05-10 


13.2 


289 


8."5 


11.91 


-0.36 ±0.30 


53 


20550+1656 


0670140101 


2011-10-28 


58.2 


161 


11. "7 


9.13 


-0.52 ± 0.04 


2900 



Note. — Separation refers to the projected nuclear separation between the interacting galaxies; the first is angular separation while the second 
is physical separation. HR is the hardness ratio. Note that exposure times are after applying filtering for cosmic rays. 
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with an abrupt spike in counts above the base threshold in order to minimize contamination from cosmic 
rays and other short, time- variant noise sources. As the EPIC MOS cameras have a smaller pixel scale than 
the EPIC PN by nearly a factor of four, we first tried creating images and radial profiles using this data. 
Unfortunately, there were not enough counts in the MOS images to determine the radial profiles. There 
were more counts in the PN detector, however, so we proceeded to use this set for our analysis. Processed 
EPIC PN images for the three sources in our sample are shown in Figure Q] The total exposure time for 
IRAS 18329+5950 was 37.2 ks, with 9.4 ks usable after background filtering. The left panel of Figure Q] 
demonstrates that the nuclear regions of the two separate galaxies in IRAS 18329+5950 are resolved in 
X-rays. The two spectra were extracted from a region of radius 13' centered around a = 18:33:37.67, 5 = 
+59:53:23.52 for 18329E and a = 18:33:33.87, 5 = +59:53: 14.85 for 18329W. The background was selected 
from a region of the same size in a nearby location on the same CCD in order to minimize contamination 
from inter-chip variance. Additionally, the middle and right panels of Figure [T] show that for both IRAS 
19354+4559 and IRAS 20550+1656 the nuclei from the individual component galaxies were not distinctly 
separable. Thus only a single spectrum could be extracted and examined for each of these sources. We 
analyzed a region of radius 11" centered on a = 19:36:56.13, 5 = +46:06:26.24 with a final exposure time 
of 13.2 ks for IRAS 19354+4559. The total exposure time pre-filtering was 33 ks. For IRAS 20550+1656, 
we extracted a region of radius 19" centered on a = 20:57:24.09, 5 = +17:07:38.94 with an initial exposure 
duration of 77 ks and final time of 58.2 ks. 

As hinted at above, one avenue for discerning the probable X-ray source type is to compare the radial 
profile of the observed flux to the PSF of the detector. A lone AGN, being a point-like source, would have a 
profile that approximately matched the PSF's, whereas a starburst or a blend of starburst with an AGN would 
be more extended. We used the XMMSAS task eradial to calculate the sources' radial profile distributions 
at an energy of 2-8 keV. We chose this energy range because if an AGN is indeed present and is not heavily 
obscured, it would likely provide a greater fractional contribution to the hard band flux than the soft band as 
compared to a starburst region. The eradial task modeled the PSF at 5 keV, the midpoint of our hard energy 
range, using the built-in ELLBETA model type. This model is based around an elliptical King profile, with 
an added Gaussian to account for spokes. As the PSFs have an arbitrary normalization, we renormalized all 




Fig. 1 . — Smoothed, false-color XMM Newton images of the three objects in our sample with the spectral 
extraction region, described in ^3j overplotted as black circles. All images are scaled logarithmically over 
the range 0.5-10keV. Left: IRAS 18329+5950 Middle: IRAS 19354+4559 Right: IRAS 20550+1656. 
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of the counts per square arcsecond of the background-subtracted data to the peak surface brightness of its 
corresponding object so as to provide an easier comparison to the source. These results are shown in Figure 
|2] Our systems appear to be at least marginally extended with respect to the PSF over our energy range, if 
not greatly extended. The exception is 18329E, which is fairly consistent with being point-like. 



3.2. Spectral Analysis 



Following extraction, we performed spectral analysis using XSPEC vl2.7.0. We modeled each sys- 
tem with multiple components in combination to account for AGN, possibly obscured, and starbursts in an 
attempt to recover the observed count and energy distributions. We were, however, limited in the complex- 
ity of our models due to low source number counts, as can be seen in Table [T] Our fewest counts came 
from IRAS 19354+4559 with 53. Both sources in IRAS 18329+5950 had around 200 counts, and IRAS 
20550+1656 had over an order of magnitude greater at 2900. 

For IRAS 18329+5950, due to the low number counts in both sources, we decided n ot to bin the data 
and used Cash statistics rather than fit based on the more commonly used reduced Rvalue (lCashlll979l) . We 
took into account the Galactic absorption, w ith hydrogen column dens ity along this line of sight of A^ al = 
4.48 x 10 20 cm- 2 Jpickev & Lockmarll99o}) . using the "tbabs" model dwilms et allboooh . 183 29E was best 



fit by a power law component of spectral index T = 1.5 ±0.5 and an additional MEKAL plasma (IMewe et al. 



1985l.[l983: Eiedahl et all l995n of temperature 0.65+q }y keV to account for the soft excess shown in Figure 
[3] The final model had C = 94.78 for 115 degrees of freedom. The overall Galactic-absorption-corrected 



X-ray (0.5-10 keV) luminosity of the galaxy using this model is found to be Lx = 7.6*2 2 x 10 40 erg s 
From Table |2l we note that this source is dominated in the hard band by the power law component, and 
most likely dominated by the plasma in the soft band, although the two are roughly equivalent in this energy 
range. There is a clear soft excess component around 0.5 keV. We try to recover this with a stronger plasma 
and a steeper power law, but the fits do not converge in this scenario. 

The system 18329W, modeled in Figure 01 was also best fit by a Galactic absorbed power law and 
MEKAL plasma with C = 92.31 for 116 degrees of freedom. Here, the spectral index was T = 1.6^ 5 and 
the temperature of the plasma was 0.69 + q 2 q keV. It should be not ed that both of these models were fixed 
at solar metallicity via the prescription in lAnders & Grevessd (|1989T ). which is used by default in XSPEC. 
Like with 18329E, Table [2] indicates that 18329W derives most of its model's hard band luminosity from 
the power law and comparable amounts from both the plasma and power law in the soft band, with a total 
Galactic-absorption-corrected X-ray luminosity of Lx = 8-7^'q x 10 40 erg s . 

With fewer than a third the counts of either source in IRAS 18329+5950, the data for IRAS 19354+4559 
were also fit without binning and using Cash statistics rather than reduced \ 2 - Even so, multicomponent 
models were not well-constrained due to the paucity of counts. As such, we tried modeling the system with a 
plasma, black body, and power law component separately, each with Galactic absorption. The best fit, seen in 
Figure[5l (with C = 36.13 for 45 degrees of freedom) is the power law, with a spectral ind ex of 2.4^ q. After 
accounting for Galactic absorption through a column density A^ al = 1.2 x 10 21 cm" 2 (IDickey & Lockman 
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Fig. 2. — Radial profile plots for the three sources at 2-8 keV as well as the models of the detector's PSF 
at 5keV at the same locations. Top: IRAS 18329+5950 West (18329W, left) and IRAS 18329+5950 East 
(18329E, right). Bottom: IRAS 19354+4559 (left) and IRAS 20550+1656 (right). 
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Table 2. Model Component Luminosities 



IRAS Name 


Component 


^soft 


^hard 


^soft + hard 






log(L band /L ) 


log(L band /L ) 


log(L band /L ) 


18329+5950E 


Full Model 


^ n , if) 1 7 

'• U1 -0.17 


6 - 98 -o:n 


i—i -i r\ _i_n 1 "7 

7 30 +ul/ 




Pn\t/pr T 

r uwci j-jdw 


6 60+ 018 

°- DU -0.12 


6 Q7+ 018 


7 ia+0.18 
'• i:s -0.13 




MEKAL Plasma 


6 79+°- 16 


r 9^+0.15 
J - z - 3 -0.20 


6 81+0.15 
b - 8 1-0.20 


18329+5950W 


Full Model 


7 04+ 013 


7 07 +0 05 


7 %+ 009 

'•- )D -0.27 




Power Law 


6 60 +018 


6 93+°- 18 

°-^- , -0.12 


7 10 +018 




MEKAL Plasma 


6 73 +0 21 


r or +0.20 
J - ZJ -0.29 


6 75+ 20 


19354+4559 


Power Law 


7 06+ 017 


7 07 +017 
'■ u/ -0.24 


7 07+0. 16 
' " , '-0.24 


20550+1656 


Full Model 


7 £0+0.05 


7 qq+0.07 


8-09^06 




Power Law 


7 oq+0.07 
'• i8 -0.07 


7 QQ+0.06 






MEKAL Plasma 


7 ^7+0.05 
'•-"-0.07 


6 20+ 06 
°- zu -0.06 


7 OQ+0.06 

' --^-o.oe 



Note. — Comparison of X-ray luminosities of various model components. Luminosi- 
ties are reported in the form log(Lband/L©). All models have been corrected for Galactic 
absorption. Also provided are the individual luminosities of each of the model compo- 
nents, except in the case of IRAS 19354+4559, which has only a single power law. Errors 
are calculated from the errors on the normalizations of each model component. 




Energy (keV) 



Fig. 3. — The best fit spectrum of the eastern source in IRAS 18329+5950 (18329E), composed of a power 
law of spec tral index 1.5, a thermal pl asma component with temperature 0.65 keV, and Galactic absorption 
taken from lDickey & Lockmanl (11 9901) . The data is rebinned for illustration purposes. The bottom panel is 
the residuals between the fit and the data in units of counts per second. 
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1990D, the combined soft (0.5-2 keV) and hard X-ray (2-10 keV) luminosity was L x = 9.0^ j x 10 40 erg s" 1 



We then calculated the hardness ratio, HR, defined as 

H-S 

HR = , (1) 

where H is the total counts in the hard band and S is the total counts in the soft band. For this system, we 
find the hardness ratio to be HR = -0.36 ± 0.30. While not conclusive, this negative hardness rati o suggests 



that th e X-ray luminosity is most likely due to a starburst population or a heavily-obscured AGN (IFiore et al. 



2000). 



Finally, IRAS 20550+1656 was more complex to model due to the presence of bumps around 2 and 
4-5 keV. Because of the larger number of counts, however, we were able to bin the raw spectrum such that 
each bin contained a minimum of 20 counts. The number of counts also permitted us to use reduced x 
fitting for this system. After being unable to fit the data well with any single component model, we tried 
several multicomponent models. We found the best fit consisted of Galactic absorptio n of a power law and 



plasm a. Again, we corrected for the Galactic column density A^ al = 1.10 x 10 21 cm 2 (|Dickey & Lockman 



1990). The spectral index of the best-fit power law was T = 1.5+q 2 . The best-fit plasma temperature was 



o.i- 

0.67 ±4 keV. Universal solar abundances could reproduce most of the data, but were not able to explain the 
bump-like features. We then fit instead with a metal variant plasma, VMEKAL in XSPEC, to determine if 
any abundance patterns could match the observed spectrum. Testing each element individually, we found 
that the bump around 2 keV could be well-fit by super-solar silicon. We might expect to see enhanced 
levels of silicon from a starbursting region. This is because starbursts are a rich environment of massive 
stars, which become Type II supernovae and can pollute the region with alpha elements. To test this idea, 
we next fit the system with a model where all of the alpha elements were tied to the silicon value. This 
provided a better reduced x 2 , 0.8 for 109 degrees of freedom, than the solar metallicity plasma as well as the 
plasma with solely elevated silicon. It should be noted that we also fit with both alphas independent from 
one another as well as all metals as free parameters, but our data was unable to provide firm constraining 
power in either of these instances. The final fit, with the alphas fixed to the silicon value, has an abundance 
log ( = [Si] = 0.4 ± 0.2, and is shown in Figured 



V Si © 

We are still left with the prevalent bumps between 4-5 keV, as well as several smaller ones. If we fit the 
larger feature as a Gaussian emission line, the line has a central value of E = 4.3 ± 0.4 keV and line width 
of (j = 0.4^ 3 keV. If this feature is an emission line from the source, it is most likely Ca XX at 3.01 85 A, 
or ~4. 1 keV. We found no r eports of this line in emission in the literature, though it may have been seen in 
absorption in a few systems (ITombesi et al. 2010). Due to the tenuous detection of this line and the marginal 



improvement to our fit after its addition, however, we report the model without this component as our best 
fit. 

Overall, the observed luminosity of these conditions between both the soft and hard X-ray bands (0.5- 
10 keV) was Lx = 4.8^q | x 10 41 erg s -1 after correcting for Galactic absorption. These parameters are 
summarized in Table [3] From Table [2j we note that this source is dominated in the hard band by the power 
law component, whereas the two components are comparable in the soft band, as with the other systems. 
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Fig. 4. — Similar to Figure[3l but for 18329+5950 (18329W). For this system, the power law index is 1.6 
and MEKAL thermal plasma temperature is 0.69 keV. 




Energy (keV) 

Fig. 5. — Similar to Figure [3] but for IRAS 19354+4559. Here the best fit model is comprised only of a 
power law, best fit by index 2.4, with Galactic absorption. 
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Fig. 6. — Similar to Figure[3]but for IRAS 20550+1656. The best-fit model here has a metal-variant MEKAL 
plasma of temperature kT ~0.67 keV with the strong bump feature just below 2 keV in the spectrum fit by 
adopting a super-solar alpha abundance tied to silicon (log ^<f^ = [Si] = 0.4). The index of the power law 
is T = 1.5. We note that the bottom panel is the residuals between the fit and the data in units of sigma, as 
chi-squared statistics were used instead of Cash statistics for this system alone. 



Table 3. Spectral Fits 



IRAS Name T Temp N% al (cm" 2 ) 

keV 10 20 cm" 2 



18329+5950 East 1.5; 

18329+5950 West 1.6; 

19354+4559 2.4; 

20550+1656 1.5; 



0.65;$;$ 4.48 
0.69^| 4.48 
12.0 

0.67!°$ n.0 



Note. — Comparison of spectral fitting parame- 
ters for the two components of IRAS 18329+5950 and 
the single sources in IRAS 19354+4559 and IRAS 
20550 +1656. All Galactic col umn densities are taken 
from lDickey & Lockm an ( 1990). 



4. Discussion 



With models in hand, we further investigate the starburst and AGN nature of our sources. 



4.1. Comparison to Ranalli et al. 2003 

We first compare our measured X-ray luminosities to those predicted by the IR/X-ray relationships 
presented in Equations (8) and (12) of R03 and reproduced as (0) and ® below, respectively: 



log 



-0.5-2 keV 



LfT 







log 



^FIR 

L 



3.( 



(2) 



log 



^2-10 keV 



= log 



/ ^FIR 



■3.62, 



(3) 



Lq J \ L 

The FIR luminosity is calculated from the FIR flux, Ffir, which in turn is estimated from the 60 and 100 
micron flux densities f u , measured in Janskys, by the following relation (He lou et al.ll 19851) : 



Ffjr = 1.26 x 10- n [2.58 



/o0/un /lOO/itm 

~Jy~ 



Jy 



■ + 



] Wm 2 . 



(4) 



Since these relationships pertain directly to starbursts and assume negligible AGN contamination, if they 
hold for our observations we can infer that our systems are most likely dominated by star formation rather 
than AGN. If AGN activity dominates, however, we would expect an excess in X-ray luminosity compared 
to these predictions. The la uncertainties on the intercepts that R03 provide for these linear relations are 
0.27 and 0.29 for the soft and hard band, respectively. 



4.1.1. IRAS 18329+5950 

For IRAS 18329+5950, the 60 and 100 micron fluxes reported in the RBGS predict L soft> R03 = 2.0+^ x 
10 41 erg s -1 and Lhard, R03 = 2.3*j'j x 10 41 erg s -1 from the above R03 linear relations. After accounting 
for Galactic absorption as in R03, our models find these values to be L so f t;mo dei = 3. 9+} 2 x 10 erg s , 
Lhard, model = 3.7+} g x 10 40 erg s" 1 for 18329E and L soft , modd = 4.2^ x 10 40 erg s" 1 , L hard , model = 4.5^ x 
10 40 erg s" 1 for 18329W. The system as a whole, then, has an unabsorbed X-ray luminosity of L so f t tot = 
8.1^3 2 x 10 40 erg s" 1 and L^d, tot = 8.2+3 Q x 10 40 erg s -1 . The upper bound of the soft X-ray luminosity is 
within lower range provided by Equation (f2]) and the hard X-rays from Equation ([3]) overpredict our results 
for the system as a whole. 
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4.1.2. IRAS 19354+4559 

R03 predicts L soft , R03 = 9.2^ ;§ x 10 40 erg s" 1 and L hard , R03 = 1 . 1+J;§ x 10 41 erg s" 1 for IRAS 19354+4559. 
Our model luminosities give L so f tj mo d e i = 4.5^ q x 10 40 er g s_1 an d Lhard, model = 4.5^ g x 10 40 erg s" 1 . As 
with the soft band of IRAS 18329+5950, both the soft and hard bands of IRAS 19354+4559 are on the lower 
end of the predictions. This should be viewed with caution, however, due to the limitations in this particular 
spectral model as discussed in 93.21 



4.1.3. IRAS 20550+1656 

Using the same method as for both systems above, we see that the RBGS-provided IR data for IRAS 
20550+1656 predict L so ft,R03 = 3.7+^ x !0 41 er g s_1 and L hardiR o 3 = 4.2^ ° x 10 41 erg s" 1 after correcting 
for Galactic absorption. Doing the same, we find luminosities of L S oft, model = l-S^'l x 10 erg s -1 and 
^hard, model = 2.9^ | x 10 41 erg s" 1 . Both of these are on the lower end of the error range of the R03 equations, 
but are still consistent. 

IRAS 20550+1656 was also the subject of an extensive multiwavelength campaign, including an X-ray 
analysis with Cha ndra that we can compare to our data from XMM Newton, as part of the GOALS project 



(llnami et al. They were able to distinguish two distinct objects in X-rays (see their Figure 5, labeled as 

"A" and "C+D") with Chandra's spatial resolution, whereas we see only a single source with XMM Newton. 
One of their two X-ray sources was claimed to be extended, covering two distinct IR bright regions seen with 
Spitzer. They report soft X-ray luminosities of La, soft = 6.6 x 10 40 erg s -1 and Lq+d, soft = 1.8 x 10 40 erg s -1 , 
which, combined, is less than half of the fit value from our model. Their hard X-ray luminosities for the two 
sources were both Lhard = 10 41 erg s -1 , which again sum to below the luminosity that we find. Using their 
Spitzer data on the individual sources, they conclude that source A is most likely a starburst, obeying the 
R03 relations, whereas for C+D they find that R03 overpredicts the X-ray values for its IR luminosity, as we 
find for our single source. Their hardness ratio for C+D is 0.05±0. 14, but this cannot be compared directly 
to ours due to the fact that the two instruments have different response functions. For the non-nuclear IR 
source (labeled "D" in their paper), their highest star formation rate (SFR) prediction is 120 M Q yr -1 . One 
possible explanation for the differences in luminosity between the Chandra and XMM observations is that 
IRAS 20550+1656 is intrinsically variable, which could be suggestive of an AGN in either source A, source 
C+D, variations in the extended emission from HMXBs, or combinations thereof. Unfortunately, with the 
lower spatial resolution of our observation we are unable to determine which of their reported sources, if 
only one, is responsible for our larger X-ray luminosity. 



4.2. Comparison to Symeonidis et al. 2011 



Symeonidis et all (1201 II) . hereafter SY11, examined a larger range in IR luminosities than R03 and 
found that the R03 relations, while functional for normal IR galaxies (i.e, those with Lir < 10 11 L Q ) tend 
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to overpredict X-ray luminosities for higher IR luminosities such as those in our sample. For higher, LIRG 



type I R luminosities, they found (using a modified version of the hard X-ray relation found in lLehmer et al. 



20101) the following relationships: 



log (^H(g). 4 , 5 , 

log (^) =log (g)_ 4 .38, 

where the IR luminosity is calculated from the IR flux. This, in turn, is estimated from the IRAS fluxes in 
the RBGS using the relation 

F m = l.8xl0- 1 kl3.48^ + 5.16^ + 2.58^ + ^^)Wm- 2 s- 1 . (7) 

Jy Jy Jy Jy 

The SY1 1 relations have a la error of 0.37 in the soft band and a 2a error of 0.68 in the hard band. 

Because our systems are only marginally within the luminosity range spanned by the galaxies used 
to create the R03 relationships, we also compare our model values to SY11. For IRAS 18329+5950, its 
RBGS IR fluxes predict X-ray components of L so f tj syh = 4.8^1 x 10 40 erg s" 1 and Lhard, syii = 7.2^ 7 x 

10 40 erg s" 1 . Our soft value is within the upper error bounds, and our hard value is consistent with this 
prediction. In the case of IRAS 20550+1656, Equations © and © give L softi syii = 9.4+g 8 4 x 10 40 erg s" 1 
and Lhard, syii = 1-4^'j x 10 erg s -1 . Our values are above both of these, but well within the errors of the 
soft and hard band relations. 

Noting that a lack of counts forced us to use an overly simplified model for IRAS 19354+4559, we 
cautiously investigate how our data on this system compares to values derived from the SYII relations. 
The flux upper limits from the RBGS predict L SOItj sYii = 7.7^° 3 x 10 40 erg s" 1 and Lhard, syii = l^of x 

10 41 erg s -1 . Our soft and hard band luminosities are lower than these predictions, but within the provided 
SY1 1 uncertainties, which is shown Figure|7] 



4.3. Summary of X-ray/IR Relations 

In our systems, the measured IR and X-ray luminosities are consistent to within 2a of the empirical 
relation in SY1 1, which is appropriate for higher-IR luminosity galaxies. The R03 relations slightly overpre- 
dict the X-ray luminosity, which may arise from the fact that they were derived using lower-IR luminosity 
systems. The discrepancy with R03 is, however, lower than 1.5<r on each individual object. We suggest that 
future work, when attempting to extrapolate an X-ray luminosity from IR values for star formation, use the 
SY1 1 relations when looking in the LIRG regime. 



Could something unique be happening in this luminosity space? Ilwasawa et al] (1201 lh analyzed the 
X-ray properties of a subsample of GOALS objects using Chandra and argued that there was a distinct 
change in the hard X-ray luminosity for objects in the range log(LiR/L Q ) = 11.0-11.73. They concluded 
that HMXBs fade in relevance as the dominating X-ray source for the especially IR-bright systems. What 
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replaces them is uncertain, but it is possible that AGN activity, which is not being traced by either R03 or 
SY1 1, becomes a noticeable contributor in this regime in some systems. 



4.4. Star Formation Rates 



With known X-ray and IR luminosities, we now investigate the star formation rates in these systems. 
As there are multiple relations available in the literature, we compare several different values. Using the 
X-ray star formation equations from R03, 



SFR(M© yr _1 ) = 2.2 x 10" 40 L . 5 _ 2 ke v and 
= 2.0xlO- 40 L 2 _iokeV, 



(8) 
(9) 



the apparent SFR of IRAS 18329+5950 as a whole, using our model luminosities, is 18 and 16 M Q yr -1 
based on the soft and hard X-ray bands, respectively. As described in CD these relations supposedly measure 
the X-ray luminosities from XRBs th at are born during episodes of star formation, thus tracing the SFR. R03 
calibrate Equations [8] and fusing the lKennicuttl 11998) IR starburst relations given by 



SFR(M yr -1 ) = 



1 



^FIR 



5.8 x 10 9 L m 



(10) 



From \T0\ we find star formation rates 42 M Q yr 1 . These values ar e in the regime o f starburst galaxies, 
though they span a range of almost a factor of three between R03 and lKennicutt Jl998b . 



Using the X-ray/SFR relations provided in Equations [8] and [9] for our X-ray luminosities for IRAS 
19354+4559, we find star formation rates of 10 M yr -1 and 9 M yr -1 , respectively. 

The observed soft (hard) X-ray luminosity of IRAS 20550+1656 gives a star formation rate of 40 (59) 
M yr -1 , while Equation ( flOl gives a value 79 M Q yr -1 . 

There is an appreciable spread in the range of these star formation rates, but the values imply that at 
least one of the ga laxies in each of the systems is a starburst. The R03 values in all of our sources are rather 
close to that from lKennicuttl (119981) . as would be expected based on the fact that R03 used them to calibrate 
their own relations. The inconsistency in these values, where the star formation rate derived from X-ray 
luminosities is less than those derived from IR luminosities, is also in the direction of what has been noted 
already. Our systems are X-ray faint for their IR luminosities in the standard R03 relations, so it is not too 
surprising that the R03 SFRs, which use their own X-ray/IR luminosity scaling, give larger rates for our 
syste ms which fell be low their best fit. The Kennicutt relations also have a systematic uncertainty of about 
30% (jKennicuttl 1 1 99 w . It is plausible that the combination of these two factors can account for some of the 
discrepancy in the two sets of rates. 

The underlying assumption in most of these relations is that all of the UV emission is from O and B 
stars, whose presence is indicative of recent star formation, and this emission i s repr ocessed by dust into 
the infrared. This is not necessarily the case, however, and iMiralles-Caballerol (|2012h observed some star 
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Fig. 7. — Comparison between the modeled X-ray luminosities of our three systems with predictions from 
IR/X-ray luminosity relations. Top: Plots of the R03 soft band (left, Equation [2]) and hard band (right 
Equation [3]> relationships. Crosses represent data from R03 best fit by the solid line and la errors with 
dot-dash lines. Stars represent our data, notably in a region where their data is much scarcer. Bottom: Plots 
of the SY11 soft (left, Equation [5]) and hard (right, Equation [6]) IR/X-ray relationships, with dot-dash lines 
being the one and two sigma errors for the soft and hard bands, respectively. Note that here our soft data are 
on the lower bounds of R03 and up per bounds for SY11, but with in error tolerances for all four cases. The 



extra data points were compiled by 



Pereira-Santaella e t al. (1201 lb . and taken from 



Miniutti et al 



Lehmeretal 



(2007) for MCG-03-34-064. lLevenson et al. 



(2010) for 
(2005) for 



NGC 23, NGC 5653, and Zw049.05 7, 

NGC 7130, and blustin et all (120031) for NGC 7469. Our systems tend to lie in a region occupied by HII and 
composite AGN/starbursts. 
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clusters with as much as 15% unobscured UV radiation. This finding implies that our values, which do not 
take into account any escaping UV flux, are only lower bounds on the total star formation rate. Given our 
large uncertainties, however, this may be a negligible contributing factor. 



4.5. AGN in LIRGs? 

We have employed several methods to determine whether or not our sample contains AGN or dual 
AGN. One indicator would be X-ray point sources. From Figure Q] we see that all three systems are detected 
in the entire XMM energy bandpass and that IRAS 18329+5950 is composed of two distinct sources. If 
these X-ray sources were due solely to the presence of an AGN, one would expect that their radial profiles 
would be similar to the detector's PSF From Figure[2]we see that, for the most part, the 2-8 keV emission is 
more extended than the model 5 keV PSF. It is impossible to say from this alone, however, whether AGN or 
other separate objects, such as LMXBs/HMXBs, are the culprits because the presence of a starburst could 
blur out any point-like AGN emission or large column depths could obscure strong sources. 

The second test was to compare the spectra of each system to starburst and AGN templates. Both 
sources in IRAS 18329+5950 and the single source in IRAS 20550+1656 were composed of a MEKAL 
plasma, typically associated with ionizing photons from OB stars and their supernovae, and a power law, 
which could b e from either an AGN or X-ray binaries. As for their spectral indices, X-ray binaries are 



expected, from lPersic & Rephaelil (|2002f) . to have power law slopes of ~ 1.2, though values in the literature 
range from 1-3. The power law components of IRAS 18329+5950 as well as IRAS 20550+1656 are on 
the upper end of this range. IRAS 19354+4559 is steeper but consistent within its error bars. If it is truly 
greater than 2, this could be due to the limitations in model complexity from low counts, the presence of a 
low luminosity AGN, or both. 

Lastly, we compared our X-ray values to predictions from two separate IR/X-ray relationships derived 
from star formation. The first set, from R03, was taken from a sample of lower-IR-luminosity galaxies 
filtered for AGN contamination. All of our systems' X-ray luminosities lie near or below the lower limit of 
what their IR fluxes from the RBGS would predict from star formation alone. We then examine the SY11 
relationships, which are argued to be more applicable for IR luminosities above 10 u Lq, the regime of our 
systems, and discover that our model luminosities are mostly well fit or on the higher end of expectations 
from star formation alone on this scale. 

In all of these tests, we do not find robust evidence for AGN activity in any of our sources, though we 
cannot definitively rule out AGN activity either. 

This, then, begs the question of why there are no dual AGN, or even single AGN, signatures in our sys- 
tems. The most likely explanation is a lack of imaging and spectral data, especially the latter. Currently, our 
systems he primarily between the R03 and SY1 1 predictions, but IRAS 18329+5950 and IRAS 19354+4559 
have large uncertainties in their luminosities. Better data would allow us to determine with more certainty 
whether or not these systems are consistent with either relation. It would also permit us to attempt to fit 
more complex models with higher constraining power. This is a severe problem with IRAS 19354+4559 
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in particular, as we could only fit it with a single power law component which clearly has a soft excess, as 
shown in Figure [5] The biggest obstacle among all our systems was perhaps the loss of large portions of 
our exposures due to flaring, particularly in IRAS 18329+5950 and IRAS 19354+4559. Additionally, the 
enhanced spatial resolut ion of Chandra mig ht be able to discern two distinct nuclear sources, as in the case 
of IRAS 20550+1656 in llnami et all J201ph 



It is also possible that any AGN present in our systems is Compton thick with large (> 10 cm ) 
column densities, hiding nuclear activity. Knowing that these system s are heavily obscur ed in the op- 
tical, we investigated the location of our systems on a BPT diagram ([Baldwin et al.lll98ll) of N[IIl/Ha 



again st 0[III]/H/3. In Figure M we show IRAS 18329+5950 and IRAS 20550+1656 (H/3: iKennicutt et al. 
2009, [NII1, [OIII1, and Ha: LMoustakas & Kennicuttl 120061) alongside data and classifications taken from 
Veilleux & Osterbrockl (|1987T) as reference, and see that they lie in the same region as the starbursts (SBs) 
and narrow emission line galaxies (NELGs, those galaxies with profiles akin to HII and LINER systems), 
similar to Figure [7] We also looked for mid-IR lines in the literature that wo uld be less subject to internal 
extinction than the optical ones. Only one was encountered (|Inami et al.ll2010|) . which suggested the system 
was a starburst rather than an AGN. 

Our data is also consistent with there being no AGN present at all. What does this imply? It could be 
that there is a problem in the merger scenario presented in §Q]on the major assembly and activation of most 
SMBHs by galaxy-galaxy interactions. Given the small number of systems reported with dual AGN and the 
larger number of attempts to find them through varied methods, there could also be an issue with this phase 
of the hierarchical growth route altogether. Some fraction, possibly substantial, of SMBHs may not accrete 
much material during mergers, in which case one might expect most of the AGN activity to occur during 
slower secular accretion processes. If that were indeed the case, then finding binary AGN in mergers would 
be quite difficult. Alternatively, perhaps intense star formation due to the merger process is able to consume 
most of the potential AGN fuel before it reaches the orbiting SMBHs and ignites them. The star formation 
processes also are quite powerful in mergers, and can possibly dwarf lower luminosity AGN activity that 
occurs as the SMBHs begin to accrete matter. Or, maybe the lifetime of dual A GN activity is extremely 
short and thus hard to observe, as some simulations (IVan Wassenhove et alJl2012T) have suggested. It could 
be that one of our assumptions in creating our sample selection criteria is incorrect as well. For instance, 
Class III mergers, contrary to our initial hypothesis, may be too early in the merger process to have AGN 
activity in the systems. Unfortunately, our data cannot discriminate between these scenarios and we are 
unable to definitively state which, if any, is the true culprit at this time. 



5. Conclusion 

We present the results from a search for dual AGN. Based on the environment arguments presented in 
§Q] we selected for nearby LIRGs (tending to have lower Lir than those observed at higher redshifts) with 
galaxy separations near the resolution limit of the XMM Newton and with evidence of interaction. The only 
system in which two distinct X-ray sources are resolved is IRAS 18329+5950. 



Table 4. IR and X-ray Luminosities 



IRAS Name 


^FIR 


^IR 


R03 so f t 


SYll so f t 


R03hard 


SYllhard 




log(L F iR/L ) 


log(Li R /L©) 


log(L so f t /L ) 


log(Lhard/Lo) 


log(L so ft/L ) 


log(L h ard/L ) 


18329+5950 


11.39 


11.65 


7.71 ±0.27 


7.10±0.37 


7.77 ±0.29 


7.27 ±0.68 


19354+4559 


11.06 


11.85 


7.38 ±0.27 


7.30 ±0.37 


7.44 ±0.29 


7.47 ±0.68 


20550+1656 


11.66 


11.94 


7.98 ±0.27 


7.39 ±0.37 


8.04 ±0.29 


7.56 ±0.68 



Note. — Predicted X-ray luminosities based o n IR and FIR lumi nosity relations. All IR luminosities are 
calculated from the fluxes presented in the RBGS (ISanders et al.ll2003l) . Here and throughout the paper, the soft 
band is 0.5-2 keV and the hard band is 2-10 k eV. Also note that the errors on SY1 lhard are 2er whereas the rest 
are la, or approximately so in the case of R03 (IRanalli et al.ll2003r) . 
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Fig. 8.— BPT diagram from lVeilleux & Osterbrockl (1198711 with IRAS 18329+5950 and IRAS 20550+1656 
added. They occupy the same parameter space as narrow emission line galaxies and starburst systems. 
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If we take the predicted starburst-originating X-ray values from lRanaili et al.1 (12003). the best fit mod- 
els for IRAS 18329+5950, IRAS 19354+4559, and IRAS 20550+1656 are all somewhat underluminous in 
X-rays if we, as they did, do not correct for intrinsic absorption. Given our systems' uncertainties, however, 
they mostly rema in within la of the pre dicted value for both the soft and hard bands. If we use the rela- 



tions presented in lSymeonidis et all (1201 if) for star-formation/X-ray instead, reproduced here as © and ©, 
however, we find that our models are well fit, or even slightly overluminous in X-rays for a star formation 
origin, though within 2a of the relations they provided. 

Under these circumstances, an optimistic interpretation of IRAS 18329+5950 would categorize it as a 
dual AGN candidate with hybrid starburst activity, since two separate objects are seen at the correct distances 
from one another to indicate heightened nuclear activity. The data for each of our systems, however, is not 
of sufficient quality to find or conclusively rule out the presence of AGN. Regardless of the dual AGN status 
of our systems, the l uminosities here may provide some additional support for the implications found in 



Iwasawa et al. 



(2011) and SY11, namely that there appears to be a characteristic change in X-ray behavior, 
and thus possibly a change in the sources of X-rays, in the regime of lower luminosity LIRGs, as well 
as a diminishing effect from star bursting as predicted by R03. These results could be improved upon by 
searching for other nearby LIRG systems with resolvable nuclear regions. If they continue, as the systems 
presented in this paper, to be deficient in expected X-ray luminosity, then we will either need to rethink 
our IR/X-ray relationships based on star formation in this regime or quantify the contribution from AGN 
contamination. 

We would like to thank the XMM Science Team for their invaluable advice with data processing and 
useful assistance. We are also extremely grateful to Scott Adams, Joe Antognini, Ben Shappee, Garrett 
Somers, and Obright Lorain for help along the way. 
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